function r = dx(arg)
  global d1x;
  r = d1x * arg;

  global dtest;

  if dtest
    rfft = fft2(r);
    rfft = rfft .* conj(rfft);
    keyboard();


  end
end
